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A first-principles investigation of the electron-phonon interaction in the recently synthesized os- 
mium dinitride (OSN2) compound predicts that the material is a superconductor. Superconductivity 
in OsN2 would originate from the stretching of covalently bonded dinitrogen units embedded in the 
transition-metal matrix, thus adding dinitrides to the class of superconductors containing covalently 
bonded light elements. The dinitrogen vibrations are strongly coupled to the electronic states at 
the Fermi level and generate narrow peaks in the Eliashberg spectral function a 2 F{u). The to- 
tal electron-phonon coupling of OsN2 is A = 0.37 and the estimated superconducting temperature 
T c m 1 K. We suggest that the superconducting temperature can be substantially increased by hole 
doping of the pristine compound and show that T c increases to 4 K with a doping concentration of 
0.25 holes/OsN 2 unit. 

PACS numbers: 74.10.+v;74.25. Jb;74.25.Kc; 



Metallic compounds containing light elements such 
as H, Li, and B have attracted considerable atten- 
tion recently due to their potential superconducting 
properties'^ Within the weak-coupling BCS theory, 
high-frequency phonons due to the presence of atoms 
with light masses ensure a large prefactor in the BCS 
formula for the superconducting critical temperature T c . 
Thus, even a moderate electron-phonon coupling can 
yield a sizable T c 4^ The discovery of superconductivity 
in MgB2fi B-doped diamond' B-dopcd silicon' and cal- 
cium and ytterbium graphite intercalated compounds 9 
confirms this picture and extends it by showing that 
strong covalent bonds between light atoms can provide a 
sizable contribution to the electronic density of states at 
the Fermi level, under appropriate "doping" conditions. 
In MgB2, the Fermi level crosses the covalent a bonds 
formed by boron atoms. Such states are partially empty 
as a consequence of the lowering of the ir bands, caused by 
the Mg 2+ attractive potential felt by the B-ir electrons'^ 
In B-doped diamond' substitutional boron atoms pro- 
vide hole doping to the C-C sp 3 covalent bonds. The 
strong C-C bonding allows the structure to remain stable 
even at high doping. In electron-doped graphite interca- 
lated compounds, the Fermi level crosses the graphitic 
C-7r band and the intercalated band. 

Nitrogen follows boron and carbon in the first row of 
the Periodic Table and is characterized, in its elemental 
form, by a strong triple bond in the low pressure molecu- 
lar phases and by covalent single bonds in the nonmolec- 
ular phase stable at pressures exceeding a megabar'i 
As a consequence, molecular phases are insulating and 
the nonmolecular phase is semiconducting'^ In analogy 
with the boron and carbon-based superconducting com- 
pounds described above, search for superconductivity in 
nitrogen-based systems requires the identification of com- 



pounds where covalent bonds between nitrogen atoms 
persist in a stable form in the presence of doping species 
and of a resulting metallic state. 

To our knowledge only OsN2 , a member of the family 
of late transition-metal nitrides synthesized recently at 
high pressure and temperature starting from their con- 
stituent elements ) 12 i 13 ' 14 i 15 fulfills the above criteria. The 
compounds have been obtained by subjecting the parent 
metal to extreme conditions of pressure and tempera- 
ture in a nitrogen embedding medium, in a diamond- 
anvil cell. Interest in these compounds has resided so far 
in their large bulk modulus, which suggests supcrhard 
mechanical properties. Ab initio calculations show that 
among the three compounds synthesized so far (PtN2, 
IrN2, and OSN2) only OSN2 has a metallic character , 16 i 17 
in agreement with the experimentally observed absence 
of first-order Raman peaks in this compound^ Cova- 
lently bonded dinitrogen (N2) units are preserved in the 
marcasitelike crystal structure of OsN 2 (Ref. [l6|) [see 
Fig. [Ha)] i which makes this compound an obvious can- 
didate to investigate superconductivity in nitrogen-rich 
systems. 

In this work, we investigate the superconducting prop- 
erties of OSN2 in the framework of a phonon mediated 
pairing mechanism. We performed ab initio calcula- 
tions of the Fermi surface, electronic bands, phonon dis- 
persions, and electron-phonon couplings in OsN2, and 
demonstrate that the high-frequency modes originating 
from the covalently bonded N2 units are strongly coupled 
to the electronic states at the Fermi level and would give 
rise to a T c of about 1 K. We also show that the super- 
conducting temperature can be greatly increased by hole 
doping. 

The calculations were performed with the QUANTUM 
ESPRESSO package^ employing density functional 
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theory and the Pcrdcw-Burkc-Ernzerhoff exchange- 
correlation functional^ An ultrasoft pseudopotential de- 
scription of the ion-electron interaction,— with Os 5s 
and bp semicore electrons included in the valence, was 
used together with a plane-wave basis set for the elec- 
tronic wave functions and the charge density, with en- 
ergy cutoffs of 40 and 480 Ry, respectively. The dynami- 
cal matrices and the electron-phonon coupling constants 
A were calculated using density functional perturbation 
theory (DFPT) in the linear response regime. 18 i 21 The 
electronic Brillouin zone (BZ) integration in the phonon 
calculation was sampled with a 20 x 16 x 30 uniform k- 
point mesh. The electron-phonon coupling was found to 
be converged with a finer grid of 26 x 22 x 40 k points and 
a Gaussian smearing of 0.006 Ry. The dynamical matrix 
was computed on a 2 3 mesh of phonon wave vectors q. 
The phonon dispersion was then obtained on a finer 8 3 
q mesh by Fourier interpolation of the real space inter- 
atomic force constants. In this way, A is calculated over 
a 8 3 q-point mesh. 




FIG. 1: (Color online) (a) Crystal structure of OsN2 
(isostructural to marcasite). The space group is Pnnm, with 
osmium atoms (gray) in the Wyckoff sites 2a and nitrogen 
atoms (blue) in the 4g sites (Ref. [Toh ■ (b) Brillouin zone 
showing the high-symmetry directions used in Figs. [2] and [3] 

In the energy range shown in Fig. [21 the total elec- 
tronic density of states (DOS) is essentially determined 
by Os 5d and N 2p orbitals (Fig.O right panel). At the 
Fermi level, the N 2p orbitals contribute with about 20% 
of the total DOS. The projection on the atomic orbitals 
also shows that 92% of the nitrogen contribution at 
is due to the N p XtV states and the remaining 8% to N 
p z states. The N p x ,y orbitals lie in the plane containing 
the N-N units, and are thus directly involved in the for- 
mation of the N-N covalent bond in OsN2. Integrating 
the DOS in a window of energy close to E^ (between Ep 
and 1 cV below Ep), we find an antibonding character 
for the electronic states on the N-N units. This is con- 
sistent with the considerable weakening of the N-N bond 
in OsN2 with respect to the molecular triple bond, and 
is confirmed by the large reduction of the N-N stretching 
frequency from 2300 cm -1 in the molecular state to 500- 
800 cm -1 in the compound. A non- negligible coupling of 
the electronic states close to E^ with the N-N stretch- 
ing vibrational modes can thus be anticipated based on 
simple band-structure considerations. It is interesting to 
remark that the presence of the Os framework is cru- 



cial to the presence of a finite nitrogen component in the 
DOS at Ei?. A band-structure calculation for a pure ni- 
trogen system obtained by removing the Os atoms from 
the crystal structure of OsN2 gives an insulating solid, in- 
dicating that the nitrogen component at E^ arises from 
the coupling of the N-N units with the transition-metal 
framework, through rehybridization and/or charge trans- 
fer. 
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FIG. 2: (Color online) Left panel: electronic bands of marca- 
site OsN2 at ambient pressure along the high-symmetry lines 
shown in Fig. [lib). Right panel: electronic density of states 
and its projection onto the Os 5d (blue) and N 2p (red) or- 
bitals. The Fermi energy is set to zero. 

Due to the presence of stiff N-N bonds, the calcu- 
lated phonons of OsN 2 , shown in Fig. [3j can be divided 
into three main groups: a low- frequency group (up to 
200 cm -1 ) involving mostly the Os sublattice, an inter- 
mediate group corresponding to the librational modes 
of the N-N units (between 250 and 600 cm -1 ), and a 
high-frequency manifold corresponding to rotation and 
stretching of the N-N units (above 600 cm -1 ). The 18 
phonon modes of the marcasite structure belong to eight 
irreducible representations. Four of these representations 
(Bi g , B2 g , B3 g , and A s ) are associated with nitrogen dis- 
placements only, with osmium atoms at rest, and are 
highlighted with symbols in Fig. [3l Representations B2 g 
and B3 S have one mode each and involve dinitrogen vi- 
brations along the z axis, while representations A g and 
Bi g have two modes and involve vibrations along the x 
and y directions. In B2 9 and A g both dinitrogen units 
oscillate in phase, while in B^g and Bi g the N-N units 
vibrate in counterphasc. As we can see from Fig. [31 the 
nitrogen light mass and the covalent N-N bond ensure a 
high frequency for these modes, with A g and Bi g modes 
ranging between 640 and 825 cm -1 . Raman peaks in 
this frequency range have been observed experimentally 
also in PtN 2 and IrN 2 and have been associated with 
the stretching of N-N units i 12 d 3 i 14 As a confirmation of 
the planar (xy) nature of bonding in N-N, we note that 
the rcy-polarized phonons in the A g and Bi 9 representa- 
tions are higher in frequency with respect to z-polarized 
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modes. 
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FIG. 3: Phonon dispersions along high-symmetry directions 
of the Brillouin zone and corresponding phonon density of 
states (right panel). 

Within DFPT,— the electron-phonon interaction for a 
phonon mode v with momentum q can be calculated as 



Aqjy — 



w*„N(E F )N k 
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where the sum is over the Brillouin zone. The matrix 



element is Sk„ k+qm = (kn\dV/6uq V \k + qm)/ A /2w qi/ , 
where u q! , is the amplitude of the displacement of the 
phonon and V is the Kohn-Sham potential. The electron- 
phonon coupling is calculated as a BZ average over the 
phonon wave vectors A = ^ q „ A q „/7V g . The Eliashbcrg 
spectral function a 2 F(uj) is defined as 



a 2 F{u) = — A qy w qt/ <5(w 



qi/ 



(1) 



and allows to compute A(cj) = 2 dto' 'a 2 F(ui') jio' . 

Figure Ufa) shows the Eliashberg spectral function 
a 2 F(ui) calculated from Eq. fl}. Three separate 
contributions to the electron-phonon interaction can 
be distinguished and attributed to the low-frequency, 
intermediate-frequency, and high-frequency phonons, re- 
spectively. Two high-frequency peaks are particularly 
strong and well resolved, and are associated with the 
contribution of the two A g phonons, while the contri- 
bution of the B 2ff and B 3g phonons accounts for most 
of the electron-phonon interaction in the intermediate- 
frequency range from 500 to 600 cm . The integral 
A(w), represented by a dashed line in Fig. BJa), shows 
that the low-frequency phonons that involve mostly os- 
mium atoms account for a contribution of A = 0.17. The 
high-frequency phonons associated with the stretching 
of the covalently bonded N-N units contribute with an 
equivalent amount, which brings the total A for OsN 2 
to 0.37. More insight about the nature of the electron- 
phonon interaction that leads to such a large contribution 
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FIG. 4: (Color online) (a) Eliashberg function a 2 F(u>) (con- 
tinuous line) and integrated coupling X(uj) (dashed line) of 
OsN2. The shaded regions are the B2 9 , B3 9 , and A 9 contri- 
butions to a 2 F. (b) Calculated Fermi surface of OsN2. The 
Fermi surface consists of three electron pockets located close 
to the zone center, only one of which is visible (green pocket), 
and of four hole pockets, two of them centered at Z and two 
at the S point, (c) Fermi surface upon distortion of the lattice 
along an A 9 phonon. 



to A from N-N bonds can be obtained by analyzing the 
changes of the Fermi surface (FS) that arise upon dis- 
torting the lattice along the relevant phonon modes. In 
particular, we concentrate on the A g modes that show the 
highest a 2 F{w) values and dominate the high-frequency 
contributions. In Fig. HJb), we compare the FS of the 
undistorted OsN2 with that of a distorted OsN2 crys- 
tal obtained by changing by 4% the distance between 
nitrogen atoms in the N-N units along the xy-polarizcd 
A g mode. The most relevant change in the FS is the 
migration of electrons from the bands that contain the 
hole pockets at Z to the bands that contain the electron 
pockets near T. Interband electron transfer is generally 
associated with a large electron-phonon coupling, which 
is consistent with our finding of a large contribution of 
the A g modes to A. 

The superconducting critical temperature can be esti- 
mated using the McMillan formula: 22 



1.2 



■ exp 



1.04(1 + A) 
A- A t*(l + 0.62A) 



(2) 



where /j,* is the screened Coulomb pseudopotential and 
(uj\ n ) = 280 K is the logarithmically averaged phonon 
frequency. We obtain 0.3<T C <1.2 for 0.13> /i*>0.08 
(with T c sa 1 K at the widely accepted value fj,* =0.1, see 

Rcfs. mmm. 
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A careful examination of the electronic DOS of OsN2 
(Fig. \2§ suggests that hole doping could further enhance 
T c . Hole doping in OsN2 would, in fact, lower the Fermi 
level toward a region of higher electronic DOS and would, 
at the same time, stiffen the N-N bonds by partially 
emptying the antibonding states below Ep. The re- 
ported synthesis of several nitrides from different tran- 
sition metals (Pt, Ir, Os, Pd) suggests that the synthesis 
of transition- metal nitride alloys, i.e., of compounds with 
N-N units inserted in a matrix of mixed metal composi- 
tion, is not impossible. Alloys with different composi- 
tion allow a tuning of the electronic DOS, as observed in 
transition-metal alloys^ 




a) (cm 1 ) 

FIG. 5: Eliashberg function a 2 F(uj) of hole- doped OsN2 (con- 
tinuous line) and undoped OsN2 (dashed line). 

In order to explore the consequences of hole doping, 
we carried out ab initio calculations of OsN2 with a hole 
doping of 0.5 holes/unit cell, corresponding, e.g., to a 
hypothetical alloy with Oso.75Reo.25N2 composition. As 



expected, the DOS at ~Ep increases about 2.4 times with 
respect to the undoped case. The electron-phonon cou- 
pling matrix elements remain approximately the same in 
the doped and undoped cases, but the phonon frequen- 
cies associated with the N-N modes increase by about 200 
cm -1 (see Fig. [5]), which confirms the strong coupling of 
these modes with the electronic states close to the Fermi 
level. The frequency increase causes an increase of {ui\ n ) 
to 310 K. The total electron-phonon coupling parameter, 
A, increases to 0.49, leading to a superconducting critical 
temperature of ~ 4K for doped OsN2. 

In conclusion, we predict that OsN2 is a supercon- 
ductor, and that its superconducting properties are con- 
nected with a strong coupling between the stretching 
modes of the covalcntly bonded N2 units with the elec- 
tronic states at the Fermi level, similar to what has been 
observed in a number of boron and carbon-based com- 
pounds, including MgB2. We predict an enhancement of 
the superconducting temperature by doping OsN2 with 
holes, which we believe can be achieved experimentally 
by synthesizing the nitride starting from a hole-doped Os 
alloy. We hope this work will stimulate the experimen- 
tal search for other members of the dinitridc family with 
metallic character and potential superconducting prop- 
erties. 
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